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ABSTRACT 

This paper presents a technique for blind multiuser detection in asynchronous direct sequence-code division 
multiple access (DS-CDMA) systems over Nakagami-m fading channels with impulsive noise. A new M- 
estimator is proposed and analyzed for robustifying the detector. A closed-from expression for average error 
rate of DPSK signals is derived. Average probability of error is computed to evaluate the performance of the 
new M-estimator based blind multiuser detector in comparison with the linear decorrelating detector, Huber 
and Hampel estimator based detectors. Simulation results show that the new M-estimator based detector 
performs well. 

Keywords: CDMA; impulsive noise; M-estimator; multiuser detection; Nakagami-m distribution; probability of error. 



I. INTRODUCTION 

Blind multiuser detection can be used to reduce the 
inter-symbol interference (ISl) and also to improve 
system performance [1]. Subspace-based blind 
multiuser detection technique is known for its better 
performance among all the other blind methods. In 
[2], an improved subspace-based robust blind 
multiuser detection technique for synchronous and 
asynchronous CDMA systems over non-Gaussian 
channels is presented. A new adaptive algorithm for 
blind multiuser detection of coherent BPSK signals in 
DS-CDMA systems operating under non-Gaussian 
impulsive noise is presented in [1]. The problem of 
robust multiuser detection in non-Gaussian channels 
has been addressed in the literature [3], which was 
developed based on the Huber and the Hampel M- 
estimators, respectively. Differential non-coherent 
data detection for CDMA over flat-fading non- 
Gaussian channels with impulsive noise is presented 
in [2], [4]. Robust multiuser detection in 
synchronous DS-CDMA system with MRC receive 
diversity over Nakagami-m fading channel is 
presented in [5] by assuming that the modulation is 



binay PSK (BPSK). But, the differential modulation 
must be used to overcome the phase ambiguity [1], 
[4]. Recently, [6] considered the robust blind 
multiuser detection for synchronous DS-CDMA 
systems over Nakagami-m fading channels. 

This paper extends the work of [6] to asynchronous 
DS-CDMA systems over Nakagami-m channels in 
the presence of heavy-tailed impulsive noise by 
employing subpace-based blind multiuser detection 
technique. An approximate expression for average 
probability of error of an M-decorrelator derived in 
[6] by using an approximation to Marcum-Q 
function. In this paper, a closed-form exact 
expression for average error rate is derived using 
[Eq. 3, 7]. This expression is used to compute the 
probability of error of linear decorrelating detector, 
the Huber and the Hampel estimator based detectors, 
and the newly (modified-Hampel) proposed M- 
estimator based detector. Simulation results shows 
that the new M-estimator based detector outperforms 
the linear decorrelating detector, Huber and Hampel 
estimator based detectors. 
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The remaining part of the paper is organized as: DS- 
CDMA system over multipath fading channel in 
impulsive noise environment is presented in Section 
II. In Section III, new M-estimator proposed to 
robustify the detector is presented. Section IV 
presents the subspace-based robust blind multiuser 
detection technique. An expression for average 
probability of error of M-decorrelator is derived in 
Section V. Section VI, presents the simulation results 
and finally, conclusions are drawn in Section VII. 

II. SYSTEM MODEL 

An asynchronous DS-CDMA system shared by L- 
users, where each user transmits information by 
modulating a signature sequence over a single-path 
Nakagami-m fading channel, is considered. The 
received signal over one symbol duration can be 
modeled, when the users are symbol and chip 
asynchronous, as [4]: 



L M-l 



rit) = mlY^^bi (i)ai (t)e^ l^siit- iT, - r, ) + «(0 (1 ) 

U=l !=0 J 

where 

denotes the real part, M is the number of data 
symbols per user in the data frame of interest, Ts is 
the symbol interval, «/(0 is the time-varying fading 

gain of the I* user's channel, djit) is the time-varying 
phase of the I* user's channel, bi{i) is the bit of the 
user, Si (t) is the normalized signaling waveform of 
the user, v is the propagation delay and n(t) is 
assumed as a zero-mean complex non-Gaussian 
noise given by [3] 



f = (l-£mO,V^) + £^(0,KV^) (2) 



with u>0, 0 < f < 1, and s: > 1 . 



component, with s representing the probability that 
impulses occur. At the receiver, the received signal 
r(f) is first filtered by a chip-matched filter and then 

sampled at the chip rate, 1 / r^. . The resulting discrete- 
time signal sample corresponding to the chip of 
the ;* symbol is given by [4] 

f2~'T,+(n+i)r 
r„(0 = j— I r(t)e-J''^'dt,n = l...N . (3) 

iT,+nT. 

Assuming that the fading process for each user varies 
at a slower rate that the phase and amplitude can 
taken to be constant over the duration of a bit, (3) 
simplifies to [4] 

rj ) = ^ [si (i)bi ii)h) m+gi {i-\)bi {i-l)h) d)] + Wj (i) 
(4) 

where 



(l-7r,)a'. 



J-m, 



j<mi 
j = mi 



(5) 



;r/fl'. , +(l-;r/ )J. , / > mi 



hj(l)D 



with nil = 



Kia 



= mi (6) 
j>mi 



- nil 



w = I n(t)e J'"''dt. 



(7) 



T T 
Let/:(i) □ [ro(0,-/Af-i(0] , D \_WQ(i\...,Wf^_x{i)\ and 



Here X(0, v ) represents the nominal background 

2 

noise and the K(0, ktu ) represents an impulsive 



^(0 0 -^[h(.i)g\(i)^-Hii)gL(.i)f ■ By Stacking k 

successive data samples, the following quantities can 
be defined [2], [4], 
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resist outliers without preprocessing the data [9]. In 
M-estimates, unknown parameters ^j,^2'- -^l 
solved by minimizing a sum of function p{-) of the 
residuals [4] 



(k+l)N 



(8) 



where 



with 



H{0) 



and 



m^) H(o) 0 

0 



H(l) H(0) 



(10) 



h (0) 



.%(0) 



.%(!) 



4(0). 



/if(l) 



(11) 



(12) 



The parameter k is called the smoothing factor and is 
chosen such that the matrix //^ has full column rank. 
III. M- ESTIMATOR 

There exist in the literature a number of different 
approaches to the robust estimation problem, and the 
M-estimation is one of the most sophisticated 
approaches to this problem [8]. An important class 
of robust estimators is M-estimators and they can 



9 = ai-gmin E {p 



+ P 



9? r„(/)-Z[A]„,^,(0 

1=1 



3 r„(/)-Z[A]„,^,(/) 

/=l 



(13) 



where /? is a symmetric positive-definite function 
with a unique minimum at zero, and is chosen to be 
less increasing than square, r^, (/) and 0/ (i) are the 
and I* elements of the vectors r(i) and 
respectively, lA]^i is the n/* element of the 
matrix A [see 4 for details], and denotes 
imaginary part. Since least squares (LS) estimate is 
very sensitive to the tail behavior of the noise 
density, linear decorrelating multiuser detector is 
also sensitive to the tail behavior of the noise 
distribution. Hence, a new M-estimator [10] is used 
to deal with heavy-tailed noise. The influence 
function of new M-estimator is given by [10], [11] 



C for\C\<a 
flsgn(^) for a<\(^\<l 

2^ 



a ^ 



(14) 



for \C\>b 



2 2 

The choice of the constants a {=kv ) and h {= 2ko ) 
depends on the robustness measures derived from 
the influence function. 
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IV. ROBUST BLIND MULTIUSER DETECTION 

In robust blind multiuser detection for asynchronous 
channels the signal subspace components of the 
covariance matrix of the signal are computed 

as [2], [4] 

72 ttT 



C ^E[ri(i)ri(if] = H,(i)WlHl + Ol 



N 



(15) 



Here, the signal subspace has dimension of L{m+1). 
Then, the robust estimation in the signal subspace 
follows as 

k . 



'■{i} = y^(r-u/{i)) (16) 



(17) 



1 2 

where /i is a step-size parameter chosen as — = v . 

The estimate of the product of the amplitude and the 
data bit 0; (/) = W;/?; (0 of the /* user is given by [4] 



L{k+\) u^.]y 

Olii)= I -^^/(0,/=l,2,...,L 



7=1 



(18) 



where h 



/ ^N(k-2) 



and finally, the user 



data bit is demodulated as 



bi{i + l) = sga<-J{ 



(19) 



V. AVERAGE PROBABILITY OF ERROR OF M- 
DECORRELATOR 

The asymptotic probability of error of non-coherent 
DPSK demodulator can be evaluated by averaging 
the conditional probability [4] 



p (m (0^,* (i - 1)} < 0 |a (i) = 0, X, y j 



= 2, 



M 



T 2 



-X, 



2at 



f \ 



1 



V J 



4a: 



h (20) 



over the joint probability density function (PDF) of 
the random variables X and y [4]. Here 

-i" 



a- □ V 



N 



R 



11 



(21) 



with 



2 1 1// (u)f(u)du 

V = 2 

[\¥'iu)f{u)du\ 



(22) 



and R is the cross-correlation matrix of the random 
infinite-length signature waveforms of the L users, 
the Marcum's Q-function and Ig{-)\s the 

modified Bessel function of the first kind with order 
zero. Here, it is assumed that the random variables x 
and y are statistically independent and Nakagami-m 
distributed with PDF [12]. By following the steps as 
in [6] and using [Eq. 3, 7], the average probability of 
error of DPSK signals can be derived as (23) and 
shown at the end of the paper. In (23), 



H 



(l/4cr? )/(4m/Q+l/cr? ) 

^1 ^ 

(OT/Q+l/4a? ) 

^1 



, 2^l(-)is the Gauss hyper- 



geometric function and r( ) is the gamma function 
[13], 



N 



+ 2SNRja{\-pi)\, 



(23) 



2 2 2 2 2 

a = (1 - s)v + SKV (with o and kv , respectively, 
as the variance of in-phase and quadrature 
components of noise samples) and 



SNRi □ 



\si(i)\ 



(24) 



VI. SIMULATION RESULTS 

In this section, the performance of M-decorrelator is 
presented by computing (23) for different values of 
Nakagami fading parameter and different order of 
diversity. It is assumed that the channel model is a 
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lightly damped second-order autoregressive (AR) 
process given by [4] with the bit rate, the pole radius, 
and the spectral peak frequencies are Tb = 10 kbps, rd = 

0. 998 and fp = 80 Hz, respectively. Further, it is also 
assumed that M = 1 in evaluating average probability 
of error, (23). 

Performance of decorrelating detector with different 
influence functions is shown in Fig. 1 and Fig. 2. In 
Fig. 1, the average probability of error versus the 
signal-to-noise ratio (SNR) corresponding to the user 
1 under perfect power control of an asynchronous 
DS-CDMA system with six users (L = 6) and a large 
processing gain, N = 127 is plotted for m = 1 and D = 

1. Similarly, in Fig. 2, average probability of error is 
plotted for m = 2 and D = 2. The simulation results 
reveal that the increase in diversity order (from 1 to 
2) improves the detector performance. Simulation 
results also reveals that the proposed M-estimator 
outperforms the linear decorrelating detector and 
minimax decorrelating detector (both with Huber 
and Hampel estimators), even in highly impulsive 
noise. Moreover, this performance gain increases as 
the SNR increases. 

VII. CONCLUDING REMARKS 

Blind multiuser detection in asynchronous DS- 
CDMA systems over Nakagami-m fading channels in 
an impulsive noise environment is presented. A new 
M-estimator (modified-Hampel) based blind 
multiuser detection technique is presented and its 
performance is analyzed by deriving a closed-form 
expression for average error rate. This expression is 
used to compute the average probability of error of 
the M-decorrelator with least-squares, Huber, 
Hampel and proposed M-estimators. Simulation 
results show that the blind multiuser detector offers 
significant performance gain over the linear 
multiuser detector and the minimax decorrelating 
detectors with Huber and Hampel M-estimators, in 
impulsive noise. Effect of fading parameter and 
diversity order on the performance of decorrelator is 
also studied. 



= 0,01. ic= I00.L = 6 




Fig. 1. Average probability of error versus SNR for 
user 1 for linear multiuser detector (Least Squares), 
minimax detector with Huber, Hampel and proposed 
M-estimator in asynchronous DS-CDMA channel 
with moderate (s = 0.01) impulsive noise. 




Fig. 2. Average probability of error versus SNR for 
user 1 for linear multiuser detector (Least Squares), 
minimax detector with Huber, Hampel and proposed 
M-estimator in asynchronous DS-CDMA channel 
with highly (s = 0.1) impulsive noise. 
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